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Stochastic mixing model with power law decay of 

variance 

By S. Fedotov, M. Ihme and H. Pitsch 


1. Motivation and objectives 

The mixing of a conserved scalar c = c(t,x), advected by a turbulent flow, remains 
a problem of both fundamental and practical interest. One of the basic characteristics 
of the mixing process is the rate at which the scalar variance <J%(t) = ((c — n) 2 ) decays 
with time. Here /r is the mean value, and the angular brackets ( • ) denote an averaging 
procedure. One of the simplest and widely used mixing models is the interaction by 
exchange with the mean (IEM) (Villermaux & Devillon 1972) or linear mean square 
estimate model (LMSE) (Dopazo & O’Brien 1974, Sabel’nikov & Gorokhovski 2001). In 
this model, the scalar relaxes toward its mean /i according to the equation 


dc 
d t 


r 


(i.i) 


In the coalescence-dispersion model (CD) (Curl 1963), mixing of two particle volumes 
is described by (1) a coalescence and (2) dispersion process during one time step At. 
After that time, both particles have the same scalar value, which is equal to the mean 
of the two values before the mixing. The pdf relaxes to a bell-shaped distribution but 
deviates from the Gaussian shape primarily in the tails of the distribution (Peters 2000). 
Both models introduce a time scale r which is commonly approximated by a turbulent 
time scale. Other mixing models are based on mapping closure, Fokkcr-Planck model, 
euclidean minimum spanning trees (EMST), or Langevin model. Reviews of these can 
be found in Peters (2000) and Pope (2000). 

There are two different laws governing the decay rate of a passive scalar: 

(i) the exponential law 

<y 2 {t) oc exp{— f/r} (1.2) 

with the characteristic time r and 

(ii) the power law 

al(t)oct~ a (1.3) 

without any characteristic time scale (Lesieur 1997). 

Most theoretical models introduce a characteristic time scale and assume implicitly or 
explicitly the exponential decay rate (1.2) which is only appropriate for stationary tur- 
bulence. 

The main purpose of this paper is to study the mixing process following the power law 
(1.3) that is typical for decaying turbulence. Experimental results show that the decay of 
the variance strongly depends on the initial ratio of the velocity and scalar length-scales 
and that there is no universal decay exponent (Durbin 1982). 

It is well known that the decay exponent a depends on the low wavenumber part of 
the scalar and velocity spectrum, E c (k,t) and E(k,t) (Lesieur 1997). One can expand 
both spectra into a Taylor series: E c (k,t) = 27rfc 2 (Co + C%k 2 + 0{k 4 )) and E(k,t) = 
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2irk 2 (Bo + B 2 k 2 + 0(k A )) (Chasnov 1994). For high Reynolds and Peclet numbers, the 
scalar variance can be considered as a function of Cq,B 0 and t alone and dimensional 
arguments lead to 

a 2 c (t) cx C 0 B~ 3/5 t " 6 / 5 (1.4) 

with a = 1.2. For Cq, B 0 equal to zero, the scaling argument results in different approx- 
imate decay laws with a equals to 6/7, 10/7 or 2 (Lesieur 1997). 

Here we present a simple stochastic mixing model based on the law of large numbers 
(LLN) (Feller 1966). The reason why the LLN is involved in our formulation of the 
mixing problem is that the random conserved scalar c = c(t,x(t)) appears to behave 
as a sample mean. It converges to the mean value p, while the variance <J 2 (t) decays 
approximately as t _1 . Since the variance of the scalar decays faster than a sample mean 
(typically a is greater than unity), we will introduce some non-linear modifications into 
the corresponding pdf-equation (see Eq. (2.28) below). The main idea is to develop a 
robust model which is independent from restrictive assumptions about the shape of the 
pdf. Here we exploit the similarity of the behavior of a scalar c to that of the sample 
mean 

1 n 

Cn = -V'Cfe, (1-5) 

where Cl, ■ • ■ , Cn is a sequence of mutually independent random variables, each having a 
mean p and standard deviation a 2 . The LLN tells us that the random sum c n tends to 
the mean value p with probability one, while the variance cr 2 n = (( c n — p) 2 ) decays as 
a 2 /n. Then, any arbitrary initial pdf p(n o, c) tends to a ^-distribution S(c— p) as n — > oo. 
In the present paper the discrete increment n can be understood as the time variable t. 

The main result of this paper is the derivation of the time-discrete non-linear integral 
equation for the pdf of c 


i >, 

p(t + l,c) = J j 


t + 1 


;P (t 


0 -1 


(t + l)(c - p) - A(1 + s)(r] - p) 


1 + t — A(l+e) \ \ + t — A(1 + e) 

dsdrj , 


( 1 . 6 ) 


for t = 1,2,3,... . Here "0(f,£) is the pdf for the exchange rate and A is the mixing 
intensity (see below). The main property of this equation is that it describes the relaxation 
from an arbitrary initial distribution to a d-function 


p{t,c) — > <f(c — p) as t — » oo (1.7) 

and the decay of the variance, a 2 (t) = (( c(t ) — p) 2 ), is of the form t~ a . The case a = 1 
corresponds to the law of large numbers. 

The remainder of this paper is organized as follows. In Section 2 we derive the integral 
equation from a stochastic difference equation describing the evolution of the pdf of a 
passive scalar in time. The stochastic difference equation introduces an exchange rate 7 „ 
which we model in a first step as a deterministic function. In a second step, we generalize 
as a stochastic variable taking fluctuations in the inhomogeneous environment into 
account. In Section 3 we solve the non-linear integral equation numerically and analyze 
the influence of the different parameters on the decay rate. The paper finishes with a 
conclusion. 
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Figure 1. Mixing problem along Lagrangian path way. 


2. Mixing model 

2.1. Problem statement 

The evolution of a passive scalar c = c(t , x) is governed by the stochastic PDE 

dc , . o 

— + v(t, x)Vc = DV 2 c, 


(2.1) 


where v(t, x ) is the random velocity field and D is the molecular (or thermal) diffusivity. 
The classical problem is to derive a closed equation for the Euler one-point probability 
density function (pdf) p = p(c;t,x ) = (S(c — c(t, a :))). A detailed discussion of this, 
still unsolved, problem can be found in Pope (2000). In the present paper we consider 
the mixing problem in the Lagrangian framework by introducing a passive scalar c(t) = 
c(t,x(t)) of a particle volume moving with the velocity v(t,x(t)) (see Fig. 1). 


2.2. Stochastic difference equation 

The equation for the scalar c(t) = c(t,x(t)) of the particle volume moving with the 
random velocity v(t,x(t)) can be approximated by the stochastic equation 


J = -7(t)(c-c(i)), (2.2) 

where the exchange rate 7 (t) and the ambient concentration c(t) are random processes. 
For the constant values of 7 and c, we obtain the IEM model (1.1) with the exponential 
decay: 

“TXT — ^ = exp{- 7 f} . (2.3) 

c(0) — c 

The key feature of the present model is that c(t) is a random process and functional of 
c(t) itself. The crucial assumption that relates Eq. (2.2) to the law of large numbers is 
that the mean value (7(f)) behaves as f -1 for large t. One can show that for a Gaussian 
scalar field c = c(t,x), the mean value (7 (t)) is proportional to the ratio x/cr 2 , where x 
is the mean scalar dissipation rate, % = 2 H((Vc) 2 ) (Pope 2000). It follows from 


der 


2 

C 


d t 


= -X 


(2.4) 
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that if ( j;?(t ) a t~ a then x decays as t~ < - 1+a ' > 1 and therefore (7 (t)) oc t _1 . The appearance 
of the power law decay can be understood if we assume that 7 = A (t 0 + f) _1 


dc A 

— = (c 

df to + 1 

The solution to this simplified equation is 


2 ). 


c(t) - C = / to \ A 

c(0)-c \t 0 + t ) 


(2.5) 


( 2 . 6 ) 


Since our aim here is to relate the mixing problem to the law of large numbers, it is more 
convenient to rewrite Eq. (2.2) as a stochastic difference equation 


C n +1 - C n = - 7 n(c n - c n+ i) , 0 < c n < 1 , for n = 1, 2, 3, . . . (2.7) 


where j n and c n are assumed to be sequences of mutually independent random variables 
with the densities 


w(n,7) 

and the first moments 


cl 

dy 


P{ln < 7} , 


(7™) 


1 

J 7^(n,7)d7 

0 


A 

1 + n' 


0(n,c) = —P{c n < c} 
dc 


1 

(c n ) = J c 4>{n,c)dc = n . 


( 2 . 8 ) 


(2.9) 


Here the mean exchange parameter (”/ n ) has been chosen in such a way that it is equal 
to A/2 at time n = 1. The parameter A can be regarded as a measure of the mixing 
intensity. Since 0 < c n < 1, 7„ obeys the inequality 0 < < 1. By introducing the 

deviations from the mean /r for the concentration of the particle and its surrounding, 
respectively, 



'U'n — Cn 5 


(2.10) 


£n = C n 


(2.11) 

we can rewrite Eq. (2.7) as 




^n+l r )n{^ j n £n+l) i 

~H<u n < 1 — M , 

n= 1,2,3,... , 

(2.12) 


where is a sequence of zero mean, independent random variables with the density 
(/(n,£) and — /i < £ < 1 — /i. 


2.3. Law of large numbers and forward Kolmogorov equation 

To illustrate the connection between the mixing problem and the law of large numbers, 
consider the case when the sequence f n is stationary, that is y;(£) is independent of n, 
and the exchange parameter 7„ is a deterministic sequence of the form 


In = 


1 

1 + n 


(2.13) 


with A = 1. If we assume u\ = £1, it follows from Eq. (2.12) that U 2 = (£1 + ^2)/2,«3 = 
(£1 + £2 + &)/3 and so on. Therefore, the solution of the equation (2.12) can be written 
as a sample mean 


u n 


1 

n 




k = 1 


(2.14) 
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which tends to the mean p as n — > oo while the variance (u„) oc n^ 1 . 

The advantage of having an equation (2.12) is that we can easily derive an equation 
for the probability density function 

p(n, u) =a -—P{u n < it} , n = 1, 2, 3, ... . (2-15) 

dit 

It follows from Eqs. (2.12) and (2.13) that u n is a discrete Markov processes and the pdf 
for this process satisfies the forward Kolmogorov equation (Feller 1966) 

p(n+l,u) = ^-^ J p ^n, ( n + 1 ^ u — <p(0 d £ , n= 1,2,3,... . (2.16) 

— M 

Here we used the inverse equation 

U n = n~ x [(1 + ri)u n+ i - £ n+ i] . (2.17) 

The main properties of the solution of the Kolmogorov equation (2.16) with the arbitrary 
initial condition p(\,u) are 

p(n,u ) — * S(u) as n — > oo (2-18) 


and 

i -m 

J u 2 p(n, u)du — » n~ x as n — > oo . (2-19) 

-M 

The asymptotic behavior of p{n , u) for large is quite universal. It can be written as 

p(n,u) oc v^exp {— n5(u)} for n^> 1 , (2.20) 

where the function S(u) depends on the particular choice of the density for the 
random sequence fk- One can show that S(u) obeys the equation (Knessl et al. 1985) 

i -m 

u— S = lnM ( — ) , Mix) = [ exp {x^} y>(^)d^ . (2-21) 

d u \dzt / J 

-n 

2.4. General case: random exchange rate and non-linear equation for the pdf 

To account for the entire spectrum of time scales we now assume that the exchange 
parameter j n is a stochastic variable. It is convenient to write 

7 n = T~r — (l + £n), (2.22) 

1 + n 

where e n is the sequence of zero mean, independent random variables. Since 0 < < 1, 

it follows from Eq. (2.22) that 

1 -I- 77 

-1 < e„ < - 1 • (2.23) 

The decay of the variance of the exchange rate, <j 2 (n), can be determined as 

^(n) = ((7n-(7u)) 2 ) = 0^ 


(2.24) 
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When both e n and are random variables, we have to specify the joint probability 
density for e and £. If we denote it by ^(n, £,£), then the forward Kolmogorov equation 
for p(n, u ) takes the form 


p(n + 1, u) 


1 l + n 1 

l-M 1 


n + 1 


-i 


1 + n — A(1 + e) 


( (n + l)u — A(1 + e)£\ 

\ 1 + n — A(1 + e) J 


&(n,£,{)d£d£ , 


(2.25) 

This equation follows from the stochastic difference equation for the Markov jump process 

u n+ 1 = u n - (1 + £ n )(u n ~ £n+i) for - /x < it n < 1 - /i (2.26) 

1 + n 

and its inverse function 


(n + l)it n +i — A(1 + e)£ n +i 
1 + n — A(1 + e) 


(2.27) 


Since is determined by the random scalar dissipation rate, it is natural to assume 
that 7 n and £„ are independent. Moreover, the density for £ n must be related to p(n,u). 
The simplest choice would be to assume p(n ,£) which gives T(n,£,£) = p(n,t;)ip(n,£) 
and the generalized non-linear integral equation for p(n, u) can be written as 


p(n+ 1 ,u) 


l-M 


--1 


n + 1 


-m - 1 


1 + n — A(1 + e) 


/ (n+ l)u - A(1 + e)A 

\ ’ l + n-A(l + e) / 


p(n , ^)r/>(n, £)dsd^ 


(2.28) 

for n = 1,2,3,.... It follows from Eq. (2.24) that the variance of the exchange rate cr^(n) 
decays as n~ 2 . In the subsequent sections we assume for simplicity that the variance 
of = (e^) is constant. 

We expect that the variance of a passive scalar <r 2 (n) = ((c n — p) 2 ) = (v%) behaves as 
n~ a for large n. Taking the analytical solution of the simplified model (2.5) into account, 
we anticipate a strong dependence of a on the mixing intensity A. 


3. Results 

The integral equation (2.28) describes the relaxation of an arbitrary initial density 
distribution of a conserved scalar c to a ^-distribution for large n. Because of the non- 
linearity in c and c, the equation needs to be solved numerically for non-trivial initial 
conditions. 

The trivial case A = 0 describes the motion of an inert particle ( e.g . dye) in its surround 
and Eq. (2.12) has the simple solution u n = U\ for all n. However, if A is close to its upper 
limit, the Lagrangian particle volume experiences strong interactions with its surrounding 
and intuitively, it can be expected that the mixing intensity A has a distinct influence on 
the decay of a 2 (see solution (2.6) to simplified equation). 

The stochastic variable £ influences the exchange rate 7„ and the effect of its dis- 
tribution on the decay exponent a is desirable to understand. Formally, we can ex- 
press a as a non-explicit available function of the form a = a{\,tp{n,£)) where A = 
A(Ai/^ciPe or Re). In the subsequent sections the behavior of Eq. (2.28) for different 
parameters and initial conditions is analyzed. 
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3.1. Numerical issues 

In order to advance the density p{n + 1 ,u) in time, Eq. (2.28) must be integrated over 
the £- and e-space which results in the overall-complexity 0(N U x x N e ) for each 
time step. Here, denotes the number of grid points in the ((-direction. This makes the 
computation rather time-consuming. However, a computer code can be easily parallelized 
along the u-coordinate. 

The initial density of u relaxes to a (5-distribution with steep gradients around u = 0. 
An adaptive grid using the equidistribution principle is employed in u- and £-space in 
order to resolve the shape of the pdf properly. In the equidistribution principle (Liseikin 
1999), an initially given number of grid points is distributed in such a way that the 
condition J“ ,+1 w(u)du = (N u — l) -1 J^ M w(«)d« is satisfied. The weight function w(u) 
is a positive function measuring the variation of the solution. 


3.2. Parameters and initial conditions 


From experiments and DNS studies it is shown that the density distribution of a passive 
scalar can be reasonably approximated by a /3-distribution. The /3-distribution has non- 
zero probability in the interval [u m i n ,u max ] and has the form 


/3(u; ( u),of ) = (M max _ dmin) 1 ° b {u-u min ) a 1 {u max -u) b 1 . 

T{a)T{b) 

It is fully described by two parameters, which are functions of (u) and of 


('ll) 'Umin 

(('ll) W m in)(Wmax (^)) 1 

CL — 

^max ^min 

^ ^max ('ll) 

2 1 

J 

((^) ^rnin) (^max (^)) -] 

0 — 

^max ^min 

2 1 
L < J 


(3.1) 


(3.2) 

(3.3) 


The /3-distribution adopts a wide range of shapes: for the maximum variance of = /x(l — p) 
and zero mean a double ^-distribution at u m in and u max is approached and for small of 
a Gaussian distribution around (u) is obtained. Throughout the following simulations, 
this distribution is used as initial condition for p(l,u). 

In order to solve Eq. (2.28) we have to provide an appropriate density distribution 
for the random fluctuation e. The time-dependent distribution of e is certainly problem- 
dependent and unknown in the present consideration. Therefore we assume for simplicity 
that e has a statistically stationary distribution. A reasonable choice would be to assume 
that e is distributed following Eq. (3.1). A particular case is obtained if has a 
(5-distribution. For the following test cases we use both distributions as prescribed prob- 
abilities. The constant parameters for all simulations of Eq. (2.28) are summarized in 
Table 1. 


3.3. Influence of initial distribution in p(l,u) 

In the first numerical experiment the dependence of the decay exponent a on the initial 
conditions in u and e is investigated. Therefore, prescribed distributions with different 
initial variance for of(l) and of are used (see Table 2 and Fig. 2). For this simulation 
A is chosen to be unity and £ is in the range —1 < e < (1 + n)/A — 1. The maximum 
possible value for e is then found for n = 1. 

The decay rate of of and the decay exponent a are evaluated and plotted in Fig. 3. 
In the left column, the four different combinations with initial variance cc„(l) = 0.20 are 
plotted, the middle column shows the decay rate for the cases with cc a (l) = 0.15. In the 
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Parameter 

Description 

Value 

nEnd 

number of time steps 

100,000 

M 

mean value 

0.5 

N u 

dimension of u 

101 

N ( 

dimension of £ 

101 

N e 

dimension of e 

101 


Table 1. Fixed parameters used for the numerical simulation. 


Run p(l,u) ip(l, e) = tp(e) 


BB1 

(3(0.20) 

(3(0.500) 

BB2 

13(0.20) 

(3(0.250) 

BB3 

(3(0.20) 

/3(0.125) 

BB4 

/3(0.15) 

d(0.500) 

BB5 

/3(0.15) 

/3(0.250) 

BB6 

/3(0.15) 

/3(0.125) 


Run p(l,u) ip(l, e) = ip(e) 


BB7 

/3(0.05) 

/3(0.500) 

BB8 

/3(0.05) 

(3(0.250) 

BB9 

/3(0.05) 

(3(0.125) 

BD1 

/3(0.20) 

8(e) 

BD2 

(3( 0.15) 

8(e) 

BD3 

/3(0.05) 

8(e) 


Table 2. Combinations of initial distributions used for the parameter study: (3(c rj?) - 
/3-distribution with prescribed variance 8(() - Dirac 5-function for (. 



Figure 2. Initial conditions for p(l,u) and ip(l,e) = tp(e)\ A = 1.0. 


right column, the evolution of p(n, u) is shown for the mono-modal initial /3-distribution 
with cr^(l) = 0.05. From the bottom row of Fig. 3 it can be seen that p(n, u) experiences 
strong dynamics over the initial interval, say for n < 10. Over this interval the variance 
does not follow the power law decay. Furthermore it can be deduced that the stochastic 
perturbation £ reduces the decay rate of a \ if a \ is close to <j‘1 rnax (see solid curves). 

The influence of the initial distribution of e on the evolution of p(n , u) and consequently 
on <j„(n) vanish for large n. This is most pronounced for the cases BB7-BB9 and BD3 
(right column). Here the decay exponents for all of these cases collapse to a single curve, 
independent of ip(s). 
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0 - 2 ( 1 ) = 0.20 all 1 ) = 0.15 0 - 2 ( 1 ) = 0.05 



l°gio(«) l°gio (”) iogio)”) 


Figure 3. Evolution of cr\ and decay exponent a for the initial /3-distribution in u. 


Run 

p{l,u) 

V’(e) 

Dynl 

[c5(/i + u) + <5(1 — p, - 

u)\/2 <5(e) 

Dyn2 

/3(0.15) 

/3(0.50) 


Table 3. Combinations used to study dynamic behavior of p(n, it) for A = 0.75. 

From this simulation it can be concluded that the initially assumed distribution for e 
has only marginal influence on the decay exponent a for large n. During the simulation 
the value for a never reaches a steady state. The final convergence rate for n = 100,000 
is |da/dn| < 1.0 x 10 -6 so that a is in the range 1.50 < a < 1.65 for all test cases used 
in the present section. 

3.4. Dynamic behavior ofp(n, u) 

In the present section we are interested in the influence of -0(e) on the dynamics of 
p(n,u). The transitional behavior of <r„ can be retarded if the mixing intensity A is 
small. This effect was already discussed previously. In the present section the mixing 
intensity A = 0.75 is used. 

In order to obtain insight into the transitional behavior we use a double ^-distribution 
and a bimodal /3-distribution as initial condition for u, p(l, u) = [6(p+u) + 5(l — p — u)\/2 
and p(l,u) = /3(0.15). For -0(e) = <5(e) a simple analytical solution to Eq. (2.28) can be 
found for the first time step n = 1. For clarity the different combinations used here are 
summarized in Table 3. At n = 1 the initial distribution for the case Dynl splits up into 
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Dyn2 


Dynl 



U 



Figure 4. Transitional behavior of p(n, u ) for different initial conditions and mixing intensity. 


Run 

A 

P( 1,«) 

ip(e) [( 

-mirn £-max] 

LAI 

1.50 

0(0.15) 

0(0.20) 

[-1.1/3] 

LA2 

1.00 

0(0.15) 

0(0.20) 

[-1.1] 

LA3 

0.75 

0(0.15) 

0(0.20) 

[-1,5/3] 

LA4 

0.50 

0(0.15) 

0(0.20) 

[-1.3] 

LA5 

0.25 

0(0.15) 

0(0.20) 

[-1.7] 


Table 4. Combinations used to study the influence of the mixing intensity A on the evolution 

of the variance. 


four peaks at positions u = A(£ — y)/2 + 7 with 7 = { — /z, 1 — p} and £ = {— //, 1 — p}. In 
the succeeding time those peaks move toward the center u = 0 whereas the tails flatten. 

An entirely different dynamic behavior is obtained for the case Dyn2. The initially 
smooth density function p(l,u) transits rapidly to a mono-modal distribution and con- 
verges to the (5-pdf. With the results of this and the previous section we can summarize 
that the prescribed distribution of e influences the transitional dynamics of p(n, u), how- 
ever its effect on the decay exponent a for large n is insignificant. 

3.5. Mixing intensity A 

From experiments and DNS studies it is well-known that the decay rate strongly depends 
on the initial ratio of the velocity and scalar length-scales, £ u and i c . The free param- 
eter A is introduced in the exchange rate (2.22) in order to link the present model to 
physical processes. It needs to be emphasized that we did not attempt to reproduce any 
experimental or numerical data, amply and in great diversity available in the literature. 
Presently, we are only interested in understanding the general behavior of Eq. (2.28). It 
was concluded previously that the presumed form of ip( £ ) has only insignificant influence 
on the decay rate. However, the dynamic behavior of p(n,u) in the transitional range 
is influenced by the stochastic variable e. In anticipation of the overall weak influence 
of ijj(s) on the decay exponent we chose arbitrarily ip(e) = 0(0.20) whereas ip(e) has 
non-zero probability only in the interval —1 < £ < 2/A — 1. The parameters used for the 
five different test cases are summarized in Table 4. 

Figure 5 shows the distinct dependence of the mixing intensity on the decay of 
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Experiments on turbulent mixing show that the variance of decreases faster with in- 
creasing ratio £ u /£c, until about £„/£ c ~ 5. For ratios greater than that value the depen- 
dence is negligible. This tendency can be resembled by changing the mixing intensity A. 
The power law exponent a is approximately 1.77 for high mixing intensity and decreases 
to about 0.62 for A = 0.25 which corresponds to an indolent mixing process. From the 
present findings it is difficult to specify the exact functional relation between A and initial 
conditions used in experiments. Therefore we express the mixing intensity formally as 
A = X(£ u /£ c ,Re or Pe). 


4. Conclusion 

In the present paper we presented a robust mixing model based on the law of large 
numbers. The exchange rate decays as and links the mixing model to the law of large 
numbers. By adding stochastic fluctuations to the exchange rate, we extend the model 
to the general case, taking the fluctuations of the surround into account. 

The sensitivity of the model is analyzed for different initial conditions and parameters. 
It was shown that the initial distribution for e and u has only insignificant influence on 
the decay rate of the variance. This is to be expected because the initial distribution of 
the scalar is almost never measured in experiments, assuming implicitly that either the 
scalar is distributed initially as a double 5-function or generally disregard the dependence 
of the decay of on the initial density distribution. 

It was shown that with increasing mixing intensity the variance of the passive scalar 
decreases faster. The connection of our model to physical mixing processes is provided 
through A, which can account for initial length-scale dependences £ u /£c and other physical 
effects. 
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